Drag reduction in pipe flow by optimal forcing 
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In most settings, from international pipelines to home water supplies, the drag caused by turbulence 
raises pumping costs many times higher than if the flow were laminar. Drag reduction has therefore 
long been an aim of high priority In order to achieve this end, any drag reduction method must 
modify the turbulent mean flow. Motivated by minimization of the input energy this requires, 
linearly optimal forcing functions are examined. It is shown that the forcing mode leading to the 
greatest response of the flow is always of m — 1 azimuthal symmetry. Little evidence is seen of the 
second peak at large m (wall modes) found in analogous optimal growth calculations, which may 
have implications for control strategies. The model's prediction of large response of the large length- 
scale modes is verified in full direct numerical simulation of turbulence (Re = 5300, Re T « 180). 
Further, drag reduction of over 12% is found for finite amplitude forcing of the largest scale mode, 
m=l. Significantly, the forcing energy required is very small, being less than 2% of that by the 
through pressure, resulting in a net energy saving of over 10%. 

Eddies in turbulent flow dramatically elevate the effective viscosity of the fluid, and thereby the friction drag at the 
wall. The ever increasing pressure for energy efficiency has heightened interest in turbulent drag reduction. Whilst 
strategies for drag reduction have been developed in principle, unfortunately few have yet been realised with a net 
energy saving. 

In wall-bounded turbulent shear flows two distinct spatial scales become relevant at large Reynolds number: the 
outer scale h (e.g the radius of the pipe, or the half- width of the channel) and the inner-scale v/u T , which depends 
on the fluid kinematic viscosity v and on the 'wall'-velocity, u T = (r w /p) 1//2 , based on the wall shear-stress r w . One 
of many potential difficulties is that most of the proposed strategies are based on the control of near- wall structures, 
which both change with the Reynolds number and are extremely small at large Reynolds numbers. Turbulent skin 
friction is known to be closely associated with near- wall structures which scale on this inner scale [l| . Understandably, 
there have been many investigations to manipulate these near- wall vortical structures, for example, opposition controls 
[1] and spanwise oscillation of entire wall for plane channel flow, and [3,[5j] for pipe flow. These methods offer ^7% 
net energy savings, although a large part of the energy budget may be required to modify the flow. Passive methods, 
requiring no external energetic input other than the driving pressure, are another promising route. Streamwise aligned 
riblets have been used to produce energy savings, also around 7%, much more cheaply than via active units [y, |7|. 
Unfortunately, each riblet must scale on the inner scale, and may require cleaning should the deposition of wax or 
dirt occur. 

As an alternative route to the manipulation of individual near- wall structures, it has been suggested that many 
could be manipulated at once by means of large-scale control. The concept was proven in a nice study of the channel 
flow 01 by imposing large streamwise vortices A + « 400. The spacing of near- wall structures is typically ^100 wall- 
units [9]. Drag reduction reaching 20% was achieved, and while the effective work done to impose the vortices was not 
documented, their amplitude was only 6% of the centre-line speed. The feasibility of this strategy has been confirmed 
by an experimental investigation [l0[ in the flat plate boundary layer where large scale vortices were forced using 
vortex generators mounted on the wall. It was preliminarily shown that forcing at large spanwise scale (A + « 1000) 
is also effective in reducing the drag. 

It is convenient that large-scale vortices can induce drag reduction, as it is known they can be efficiently forced 
efficiently in laminar flows [11, 12]. It has also been realised, via the study of optimal perturbations, that the emergence 
of large scale motions is important in turbulent flows pjj EH Gil EH • Optimal perturbations consist of streamwise 
uniform vortices that optimally extract energy from the mean flow while inducing the growth of the streaks. As the 
large-scale optimal vortices scale on the outer scale, they are apt to implementation for control, as their shape is 
much more dependent on the domain geometry than on the Reynolds number. Further, the energy input necessary to 
obtain a desired amplitude of streak decreases when the Reynolds number is increased. In other words, the efficiency 
of the forcing of the large scale streaks increases with the Reynolds number. 

The scope of the present investigation is twofold. First we compute the optimal coherent streamwise vortices and 
streaks for turbulent pipe flow, in particular to determine the associated optimal energy amplification. In doing 



* Electronic address: willis@ladhyx.polytechnique.fr 
t Electronic address: yongyun@ladhyx.polytechnique.fr 
■I- Electronic address: carlo@ladhyx.polytechnique.fr 



2 



so, large growth is demonstrated, and the optimal azimuthal and streamwise wavenumbers are identified. Second, 
we ascertain that forcing at finite amplitude can induce competitive drag reduction in the turbulent pipe. In our 
simulations the circumference corresponds to A^ =1 « 1130, and the forced rolls are significatly larger than the typical 
spacing of near-wall structures. As the structure of the forced large-scale motion does not depend strongly on the 
flow rate, there is promise that a method of forcing to produce drag reduction should be practical over large ranges 
of Reynolds numbers. 

Consider the turbulent incompressible flow of a viscous fluid of kinematic viscosity v in a circular pipe of radius R. 
The bulk velocity Ub is assumed constant. Velocities are normalised by 2Ub and lengths by R. The Reynolds number 
is Re = 2 UbRjV' Following the approach used in [H, El, (3 , in order to compute optimal coherent perturbations, 
we consider small perturbations p and u, to the pressure and turbulent mean flow U = U(y)z, where y = 1 — r. The 
perturbations satisfy continuity and the linearised momentum equation: 
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where the total total normalised effective viscosity is vt(v) = 1 + E(y) with the E(y) being the eddy viscosity. We 
use for E(y) an expression originally suggested for pipe flow by Cess [17[, later used for channel flows by Reynolds 
and Tiederman Hal: 
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with R = Re / 2, B = 2 B and the fitting parameters A + = 27 and k = 0.42 chosen to better match the observations, 
more recently studied in [l9j]. The mean streamwise velocity U(y) in equilibrium with this eddy viscosity is easily 
retrieved from the mean averaged momentum equation equation 



° = - a ' p+ ib(; 
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The symmetry of the base flow allows one to consider separately perturbations of the form u(r, 0,z;t) = 
u(r,rn,a]t)e 1 ( az+rne \ where a and m are the streamwise and azimuthal wavenumbers respectively. The eigenval- 
ues and eigenvectors are found directly from the linearised primitive variable system with explicit solenoidal condition 
[26j . For all calculations considered U(y) is found to be linearly stable, in accordance with similar results obtained in 
the plane channel and boundary layer flows [l3|, [13, EH, [HI, G3 • 

The optimal response to forcing is calculated using methods described in [20[. For a harmonic forcing of the form 
/(r, m, a; t) = /(r, m, a) e 1 Qft and response u(r, m, a; t) = u(r, m, a) e 1 Qft the optimal response is given by 
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Figure [T] shows our calculations for the optimal response to forcing for several Re [27j . The optimal response i? max 
generally increases with increasing Re. As the spanwise length-scale is decreased (m increased) the response rapidly 
decreases. The largest possible response occurs for the mode m = 1. Huge response to forcing of large scale modes is 
possible. The relative difficulty of forcing motion on small scales, on the other hand, implies that considerable effort 
is required to locally manipulate structures in the neighbourhood of the wall. 

The analysis above is linear and uses the eddy viscosity assumption. A first step towards application is to verify 
that large responses are obtained in fully nonlinear turbulent numerical simulation. The pipe flow code described in 
[2l| has been used for full simulation of turbulence subject to the optimal forcing /. The numerical method uses a 
Fourier decomposition in 6 and z and finite differences in r (28[. For the following calculations the computational 
domain is of length L = 15i?, chosen to include the very large-scale motions reported in J22ll f or L up to 14. We 
choose to perform the simulations enforcing a fixed flux, consistent with previous analyses [8|, [2l| . 

Figure [2] shows time series for the skin- friction coefficient, Cf = 2(u T /U} ) ) 2 : , for three levels of forcing of the m = 1 
mode at Re = 5300. The horizontal line is the time average Cf for zero forcing. The solid line for a higher resolution 
calculation verifies no loss of the observed drag reduction for the case shown. 

The response to forcing in simulations is measured by R r = ||^||/||/||, where the averaged response field is 
u(r,9) = (u) Zt t — U(r)z and U(r) = (u)e, z ,t is the mean flow in the absence of forcing. The norm used here is 
||a|| 2 = J a • adV and ( • ) s indicates averaging over the subscripted variables. Time averages are taken over 



1400^/(2/75) time units. Measures of streak and roll components of the flow are A z 
A h = max|(u r ,U0,O)|. 
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FIG. 1: Response R m &x(a = 0, ra) to forcing, Re — 5300, 10 4 , 10 5 , 10 6 , the arrow indicating increasing Re. Inset: Re — 5300, 
m=l, response to forcing of frequency O/. The turbulent flow is most responsive to axially- independent steady forcing on the 
largest scale. 
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FIG. 2: Relative drag for several ||/|| at Re — 5300. The straight horizontal line is the unforced average. The solid line is a 
higher resolution check for the best case. 

In table [Jit is seen that the response to forcing may indeed be large, in agreement with the model. The prediction 
from the linear model is R max = 276 for m = 1, and there remains a factor 2 difference, however, with the simulations 
for the two smallest ||/||. As the time-averaged mean flow for the unforced case, obtained from numerical simulation, 
is close to that given by ([2]) and (|3]), the difference must originate from the isotropic eddy viscosity assumption. 
While at first both A z and Ah increase linearly with ||/||, streaks saturate near 19% of the centre-line speed. Streak 
saturation occurs when the roll amplitude is approximately 7% of the centre-line speed, and coincides with when the 
maximum drag reduction is achieved. 

Also in tableUis the percentage drag reduction DR = (100/C°)(C° — Cf), where superscripts indicate the unforced 
case. The power to drive the flow per unit length is P = §C/ in units p(2Ub) 3 R. The relative change in Cf for 
increased forcing is seen in figure[3j The drag reduction is up to 12.8% and the greatest net energy saving up to 11.1%, 
where E net = (100/P°)[P° — (P + W)]. For the greatest E net , the power consumption by forcing, W = \ J u ■ f dV, 
is only 1.4% of the power required to maintain the flux of the unforced flow, P°, or 1.6% of that to drive the resulting 
flow. 

The effect of forcing on the flow is shown in figure 0J Small scale axial vortices are weakened throughout the domain, 
and are almost suppressed in one half. Being largely responsible for the turbulent skin friction [23[ , this indicates that 
near- wall turbulent production is suppressed by forcing large-scale outer motions, consistent with The large scale 
forcing structure is reproduced well in the roll structure of the time averaged flow. The broad slow streak, however, 
is clearly much narrower in the nonlinear response. This narrower slow streak means that several may be packed into 
the domain (e.g. m = 2,4) before they affect each other (m = 12). The case m = 12 corresponds to the natural 
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TABLE I: Time averaged statistics for turbulent pipe flow subject to forcing. P° =3.689e-3, U c i = U(r — 0) = 0.6556. Largest 
drag reduction and net energy reduction in bold. 




mean spacing of near wall structures A + « 100. For this case the response R' is significantly reduced, and drag only 
increases as the forcing is increased. 

In summary, the model successfully predicts the large response to forcing observed in full simulation of turbulence. 
The largest response to forcing occurs for large-scale modes, which may have implications for control. Drag reduction 
of over 12% has been shown in full simulations when forcing the large-scale mode. Significantly, this was shown to be 
possible using very small input energy, being less than 2% of the energy used to drive the flow. 

Several of the alternative methods described above predict much larger drag reduction, but at much greater energetic 
expense. This is possibly linked to the small response to forcing of near- wall vortices. The net energy saving is similar 
to that here, but relies on both the large drag reduction being realised at the same time as efficient actuation. For 
useful implementation of the method proposed in this study, the induction of rolls need not be efficient. For example, 
efficiency of 50%, still using only 4% of the driving energy, might be acceptable. Note also that the response is expected 
to increase with increasing i?e, thereby reducing this cost. In practise, induction of large scale rolls is possible via 
passive actuators [24], HH or active jets [TO]. Although neither would reproduce the precise details of forcing used 
here, linear optimals are known to be robust to perturbations of the system, and responses close to those seen here 
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FIG. 4: (a) Snapshot of axial vorticity with no forcing; (b) Snapshot on same scale as (a) for 
the optimal forcing and contours of the optimal response; (d)-(f) Time averaged velocity for ||/ 
m — 1,4, 12 respectively, relative to U(r) without forcing. 
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can be expected. Slightly different responses, however, may also be beneficial — whilst the forcing in this numerical 
experiment is optimal for modifying the flow, it was not chosen to be optimal for drag reduction in the nonlinear 
regime. Thus there is scope for further enhancement. 
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Calculations were performed using a Chebyshev collocation method with N points on r £ [0, 1] . The primitive variable 
eigensystem avoids high order derivatives but has infinite eigenvalues. The eigensolver, however, returns a predictable 
number of these which are easily filtered. To be safe, 95% of the (2N — 3) eigenfunctions returned were kept for the analysis 
of optimals. At the highest Re = 10 6 with N = 250 the power spectral drop-off of the optimal mode was of 8 orders of 
magnitude 

[27] For Re T = u T R/v where u T = (ud r U(r))i \ r=R , Re = 5300, 10 4 , 10 5 , 10 6 -> Re T = 180,317,2380, 19200. 
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[28] Minimum and maximum radial spacings of points are 0.11 and 4.4 wall units, v/u T . Points are separated by 5.9 in and 
9.4 in z at r = R and L = 15 R. These values are 0.08, 4.0, 3.9, 6.3, 20 respectively for the larger calculation. 



